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ABSTRACT: The well-known practice of using the mirror-image, 
in the P,u plane, of the shock Hugoniot curve about a vertical 
line through the shocked state Phy Uy is shown to give a complete 


thermodynamic description of metals when the U(shock velocity) 
VS. Uj, relationship and o = V-1(av/aT)» 6 are known. Use of the 


experimental relationship U = Cy ot ay {a and Cy constants) and 


«a = constant, leads to a thermodynamic description which results 
in the metal appearing less compressible than if described by 
the Mie-Griineisen equation of state. Furthermore, the existance 
of an anomalous behavior of c, in the low pressure neighborhood 
(<10 to <50 kbar depending onPthe metal) of the initial state 
rules out the simultaneous existance of a Hugoniot satisfying 
the linear U vs. Wy relationship, of isentropes satisfying the 


mirror-image assumption, and of a constant value of a in this 
neighborhood. Thermodynamic functions for 16 metals are calcu- 
lated up to 2 megabars and compared with the results obtained 
from the Mie-Griineisen equation of state. 
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I. INTRODUCTION 

A probiem which is encountered quite often in 
experimental shock dynamics js the determination of the 
pressure and particle velocity behind the transmitted 
shock in a specimen after the initial shock, having 
traveled through a metal standard, has interacted with 
the interface between the specimen and the standard. 
An imoortant application of this problem is the deter- 
mination of the experimental shock Hugoniot of the 
specimen as a first step toward finding its equation 
of state, Consider a shock traveling to the right 
through a metal standard whose initial conditions (given 
by the subscript o) are u 


=O, E Vo» Py, where 


° °? 
u,E,V,P are respectively, the particle velocity, specific 
internal energy, specific volume, and pressure. A 
pressure Po of at most a few bars is considered negligible 
in present considerations since the pressures of interest 
are in the kilobar to megabar range. When this shock 
interacts with the interface a shock is transmitted into 
the specimen and a wave traveling to the left is reflected 
into the standard. The loci of P,u points that can be 
reached by the initial shock and the reflected wave in 


the standard are shown in Fig. 1. For pressures greater 


and particle velocities smaller than that corresponding 
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to the initial shock, the reflected wave is a 
reflected shock. For smaller pressures and greater 
particle velocities, the reflected wave is a reflected 
rarefaction wave. The shock equation | 

oh ae aa (1) 
relates the shock pressure Puy shock velocity U, shock 
particle velocity u,, and initial density py (=1/V,)- 
From now on the subscript H will always refer to values 
on the shock Hugoniot. Hence the measurement of p, and 
U (of the transmitted shock) in the specimen implies 


that P and u, (in the specimen) lie on the straight line ; 


n ( 
of known slope p,U as given by Eq. (1). If this is 
depicted by the dashed line in Fig. 1, then the inter- 
section with the rarefaction wave in the standard is the 
desired pressure Py and particle velocity Us behind the 
transmitted shock in the specimen (since P and wu are 
continuous at the interface). 

It will always be assumed that the effects of 
material rigidity can be neglected since the shock 
pressures are far greater than the yield point of the 
metals involved. Hence the metal can be treated as a 
fluid with an ordinary equation of state, 

The calculation of the rarefaction waves (isentropes) 


entails a knowledge of the E,P,V equation of state, 
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1 computed the temperatures along 


Walsh and Christian 
the Hugoniot and the isentropes for Al, Cu and Zn from 
the Hugoniot curves combined with the assumption that 


the specific heat at constant volume c 
2 


y and (aP/aT )., are 


constants. Walsh and Rice” noted that at least for Al 

the changes in the P,u isentropes are insensitive to 
changes in the value of (aP/2E),. Walsh, et al,” employed 
the Mie-Griineisen equation of state for which the volume 
dependence of the Gruneisen ratio was determined using 

the Dugdale-MacDonald relation. These considerations 
allowed them to compute the isentropes in the neighborhood 
of the experimental Hugoniot curve. Additional zero- 
pressure data then permitted the calculation of the 
remaining thermodynamic data of interest, and numerical 
results were tabulated for 27 metals. McQueen and Marsh? 
applied the previous theory to more extensive experimental 
data in the one-to-two megabar region for 16 metals. 

Their experimental shock measurements actually covered 

19 metallic elements. Al'tshuler, et al,2?© nave reported 
investigations of iron to shock pressures of about 5 mega- 
bars, and 3% to 5 megabars for Ag,Au,Bi,Cd,Cu,Pb,Sn, and Zn, 
The Mie-Gruneisen equation of state, in which the Gruneisen 


ratio was considered to be a constant, was used to describe 


iron. Further work! on Al,Cu, and Pb took into account the 
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electronic components of the energy and pressure. 
Consider for a moment that the specimen mentioned 
in the problem above is replaced by a vacuum (approxi- 
mating air at several bars pressure). When the plane 
shock in the metal which connects the state Po w O, 
1 O and Pips Ugy arrives at a free surface, the pressure 
in the shock is reduced to zero by a rarefaction wave 
from the surface. Let Upg be the free surface particle 


velocity. It was first postulated by Goranson® about 


1945 that for most metals 


Thus experimental values of U and Ups determine the 


Hugoniot. In fact this approximation is the first guess 


in a rapidly convergent iterative procedure for bounding? 


2,3,4 the Hugeniot when additional thermo- 


mn 


and determining 
dynamic information is given. Tabulated results?» show 
that this approximation is fairly accurate for many metals. 
For metals a further assumption (which includes 
Eq. (2) as a special case) is frequently made. The isen- 
tropic expansion from Piso Vay to 0,2u,, is approximated by 
taking the mirror-image of the Hugoniot curve about a 
vertical line through P,,,u, in the P,u plane (see, for 
example, Refs. 9 and 10). This makes it particularly ’ 


simple to calculate Py,Uy- One need have only the Hugoniot 
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P,u curve for the metal, reflect it about the appropriate 
point corresponding to the measured shock velocity in 
the metal, and find the intersection of this reflected 
curve with the straight line given by Eq. (1) (dashed 
line of Fig. 1). Under this assumption the curve of an 
isentropic compression from the point uu = 0,E,,P (R05, 
coincides with the Hugoniot in the P,u plane, but inter- 
sects it only once in the P,V plane. Hence the Hugoniot 
is not an isentrope as one might be apt to think at first 
glance, 

Since the mirror-image assumption is often used it 
is worth examining its consequences, and in the remainder 
of this paper it will be shown that this assumption and 
the experimental U, Uy relation are enough to determine a 
unique E,P,V equation of state in the metal. Furthermore 
an experimental thermal coefficient of volume expansion 
along the isobar Pi determines the temperature T, entropy 
S, and other thermodynamic properties as functions of any 
two thermodynamic variables, e.g., P and V. In Section IT 
the Hugoniot relations valid for many metals are given, 
An exact though complicated E,P,V equation of state which 
is a direct consequence of the mirror-image assumption is 
derived in Section III. Section IV contains the derivation 


of the temperature, entropy, and specific heat at constant 
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pressure. The anomalous behavior of the latter in the 
low pressure region where the Hugoniot is not well 
defined is discussed. Finally, in Section V, tabulated 
results for various thermodynamic quantities for 16 ‘ 
metals based on the mirror-image approximation are com- 
pared with results obtained by McQueen and Marsh’ who 
utilized the more sophisticated Mie-Griineisen equation 
of state. 
II, THE HUGONIOT CURVE 

Many investigators have found that for many metals 

the shock velocity is linear in the shock particle 


velocity, i1.¢e., 


U=c, + a Uy. (3) 
Theoretically, cy should be the adiabatic sound speed at 
the initial state HB Pgs Vos Uy = 0. Least-square straight- 
4 


line fits of the experimental U, Wy data’ give as the 


intercepts ec. values that are quite close to the isentropic 


° 
sound speeds at Poe Vos 


a 1 
ss 2 
0, Pav, (-32) * Utne Ue Lim vel y 2] (4) 
~ a’ 3 
OV/g Ugo PEP, YOVE 
4 Vyovy 


which were computed from Bridgeman's data by applying a 
small correction for the difference between isothermal and 
isentropic first derivatives. 


Elimination of U between Eqs. (1) and (3) yields 


6 
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Py = Polly (auyte,). (5) 
From the other two relations valid on the Hugoniot, 

vs = Py (Vy-Vy)» (6) 

EB, - EL = 2 Py (Vy-Vq)- (7) 
Viz and Ey can be written as functions of Wy? 

Vo 7 Vy = Vy Uy / (ary; + Cy), (8) 

Ey - Ey =} wa. (9) 


From Eqs. (6) and (8), the equation of the Hugoniot in 


the P,V plane is 


2 2, 
Py = 5 (¥o- ny [(2-2)¥0-2¥ | ; (10) 
the slope of the Hugoniot is 
2 
aP,, / Wy = =P, (2auyte,) (augte,)°/e <0. (11) 


If Eq. (3) holds for all shock strengths, then letting 
P,~oo in Eq. (10) leads to 


Pmax/Po = Vo/Vinin oo a/(a-1), a>l, (12) 


H 


where Prax (= 1/V.in) is the maximum deisity that can be 
obtained with a single shock starting at Po: 


III. THE E,P,V EQUATION OF STATE 
If Py = h(u,) refers to the Hugoniot pressure- 
particle velocity curve, and P(u;u,) refers to the 
| pressure-particle velocity curve for an isentropic 
expansion from a point Piys Us then the assumption that 
the expansion from the shocked state Pry Uy is the mirror 


t 
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image of the Hugoniot about a vertical line through 
Prys Uy is expressed by 
P(u;u,) = h(2uy - u), (13) 

where Uy 4s now a continuous parameter but constant for 
each isentrope. Actually Uy is a single-valued function 
of the entropy; the explicit relationship is derived in 
Section IV. By Eq. (5), 

P(usuy) = Po (2u,; - u) [a(2u, -u) + co, |. (14) 


The solution for u as a function of P is 


2u,+ %o-Uu - + ah (15) 
nH” Ba Poa (es) } 


where the positive root is taken foru < euy- Since 
aP(usu,,)/du = Pp (2au - 4 auy - Co). (16) 
the equation, 
(du) + aPav = 0, (17) 
valid on an isentrope, can now be integrated between t.:> 


points Uys Vy and u,V to give 


2p ,a(Vy-V) = ln [2 + 2a(u,,-u)/(2au, + c,) |. (18) 
Elimination of u between Eqs. (15) and (18) leads to the 


isentropic equation 


~4ap.(V-V}) 


e 2 
laP(Vsuy) = (laP, + p40,” - Poo, (19) 
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Define a new parameter 
Z = a/c, = (U - e5)/cy- (20) 
With the use of Eqs. (5), (8), and (20), Eq. (19) can 
be rewritten as 


tap (V-V,) 


tll 


ur (P,V) Cutev,P/e,2)e : 
= (2241)% en t2/(242) = wz), (21) 


where Y/(P,V) and w(z) are defined as, respectively, the 
left-and right-hand sides of Eq, (21). Using Eq. (19), 
the adiabatic sound speed c can be obtained from 

2 

p®c® = -(aP/av) 5 = Hap P + p,%c.® , (22) 
which defines the slopes of the P,V isentropes? 


The energy along the isentrope is found from 


- 
E(V3u,;) - Ey = - | P(V'5u,) av! (23) 
Va 
to be, 
2 - 2 
tal E(V5u,) - E, | - Poe (V-V55)-Vo| P(Vsu,) -P,] = 0 
(24) 


where Eq. (19) has been used. In more convenient form, 


Eq. (24) can be written as 


2 
2 
¢ = Polo PVo 
(72) O (E,P,V) = E-EQ- “i> (v-v,) 2a = 
2 
= Po%o Pav 
eo ta OW 6) ae (25) 
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Since the right-hand side of Eq. (25) is just a 


function of Uy» 


G(E,P,V) = 2°/(zt1) 2 0. (26) 
If z can be eliminated between Eqs. (21) and (26), the 
resulting equation would define the desired E,P,V equation 
of state. This may be accomplished by solving Eq. (26), 


zg? @z-O=0, 
for the single positive root z (6) and then substituting 
this into Eq. (21). In order to find the root the follow- 
ing two cases must be considered: 


(4) rf 46 < 27, then the positive root is given by 
220) (10-40 }40-0- 459}. cen 
(11) If 46 = 27, then the positive root is given by 


1 
z= (6/3 cos(.-0), (28) 
where 
a T 

0<cos0=2(27/6) s1,050<4. (29) 
Hence the isentrope characterized by z = 3 (41.e., the 
isentrope that intersects the Hugoniot when U = 4 c,) 
divides the P,V plane into two regions. Eq. (27) holds 
above this isentrope, Eq. (28) below. Therefore Eqs. 
(21), (27), (28), and (29) define the E,P,V equation of 


state, i.e., 


U(e,v) = w2| G(z,P,v) |). (30) 


10 
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IV. CALCULATION OF T,S, AND Cy 


The derivation of the thermal equation of state 
which gives the temperature as a function of pressure 


and volume proceeds in the following manner. Along an 


isentrope, 
TAS = c aT + T(aP/aT)., av = 0, (31) 
where C, is the specific heat at constant volume. Since 
c,,(8P/2E),, = (8P/aT).7, (32) 
the integral of Eq. (31) is obtained from 
T(P3z) im 
or oa ee : (33) 


The lower limit of the first integral, T(0,V), 1s found 
from the experimental coefficient of volume expansion 


along the isobar PaO, Let T(0,V) = Tof(V/V,), where T, 
is the initial temperature corresponding to EQ2Vo2Pos and 


f is an arbitrary function of WV, . From Eq. (21), 


—_ V 
g(z) ay, =] + 5 in (2z+1) - Sei) on PaO. (34) 
Define F(z) = n(0,0)/x. = f(g(z)]. Differentiation of 
(30) with respect to P gives 
(35/aP)y = | (ag /az) (az/aw)2 /aP - a /ar |az/2b, (35) 


and evaluation of the derivatives leads to 


(36) 


QB, - Eb geen). 


11 
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Substitution of Eqs. (22) and (36) into Eq. (33) gives 
the result, 


T(P32)/T, = F(z){1+aP/| p,e,(z+1)* |} ; (37) 
on the isentrope characterized by the value of z. Eqs. 
(21) and (37) determine the P,V,T equation of state through 
the parameter Zz. 

By Eqs. (25) and (34), 

ha®[E(z)-E.]/e,* = z(z-1) + 4 In(2z41) on PaO. (38) 
On the isobar P20, the entropy S can be calculated from 


T,(S-S,) = i Fay a oe) dZ. (39) 
Hence the evaluation of Eq. (39) relates S to z, namely, 
2] 2 
z Z Z aZ 
—s- [S(z)-8,] = l (2Z4+1)F(Z) * (40) 
fr) 


The positiveness of the integrand implies that S is a 
single-valued function of z. Of course this relation holds 
not only on PxO but on all paths connecting the isentropes 
Sy and 8. 

A complete E,P,V,T,S thermodynamic description of the 
metal has now been obtained. For example, given the point 
P,,V,, the parameter z, can be obtained from Eq. (21) by 
simple numerical methods. Then E,,T,, and Sy are computed 
respectively from Eqs. (25), (26), (37), and (40); in 


ce 
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general Sy necessitates a numerical integration. The 
equations of state derived above hold in the region of 
the P,V plane that lies to the right of the isentrope 


through the point Py V Ey in the direction of increas- 


o? 
ing entropy. If a Hugoniot of only finite length is 
considered (because the experimental data is limited), 
then the region of definition of the above equations of 
state is limited to the region bounded by the isentropes 
that pass through the lower and upper points of the 
Hugoniot (see Fig. 2). 


The specific heat at constant pressure, 


C= T(aS/sT), = T(dS/dz)(0z/eT),, (41) 


found by differentiation of Eqs. (37) and (40), is 


an ee, 2° (2+1)| aP+p geo” (z+1)* | 
oe pee ee 2 2 
Cy (2z+1)4 (z+1) aPtp.e, (z+1)° |dF/dz - 2aPF(z)} 
(42) 
in particular on the isobar P.O, 
(a/e,)°,0,(0,2) = 2°/| (2241) aR/az |. (43) 


Thus far the temperature variation with volume for 
PO has been given as an arbitrary function of the entropy, 
F(z(S)]. In order to find the explicit form the thermal 
coefficient of volume expansion q@ will be considered as 
constant on Pa0,Le., 


V"(aV/aT) po = a= constant (44) 


13 
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where the constant is experimentally determined. 
Integration of Eq. (44) and use of Eq. (34) leads to the 
result, 


1 2 Zz 
F(z) = 1 +a in [2 5 or 1n(2z+1) - stety | : (45) 


frei which it follows that 


Z+1 2 yA 


5 c (P,z) = 
ac, BP 1-2a oD, (z+1) (2z+1)F(z)P/ azP+p.c, z*(z+1)* |? 
(46) 
¢, (0,2) = ac®(z+1)*g(z)/a. (47) 
Since, with the use of Eq. (21), 
of F(z) = af, - in(4a) + In{4a + In Y(P,V)I, (48) 


Eq. (37) can be solved for z with the aid of hq. (48), yielding 


3 -3 
ous (ou L o 7 an(4a)F or a-in ; ] - 1. (49) 


Substitution of Eq. (49) into Eq. (21) leads to an exact 


though quite complicated P,V,T equation of state. 

Close examination of Eq. (46) shows that c, exhibits 
certain very undesirable behavior. First ¢, (P, 0) = O for 
P>O, 1.0.5 cy vanishes on the isentrope passing through the 
point E,,P,,V). Furthermore, the denominator of Eq. (46) 


vanishes on the curve C (dashed curve in Fig. 2) defined by 


P(z)=p 062° (z4+1)°/L2a“ot (z+) (22+1)F(z)e(z)-az"), (50) 


14 
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implying that c.~ oo or -oo in the neighborhood of C 


Pp 
depending on the direction of approach. This curve inter- 
sects the Hugoniot at the values of z given by the roots 
of 


z= 2a0T,(z+1)F(z)e(z). (51) 


The lower point of intersection Zs can be approximated by 


Zz, © 2aaT, <<l3 (52) 
the upper point of intersection occurs at pressures which 
are much greater than the upper limit of the experimental 
shock pressures and therefore need not be considered. 

The shock pressure at which C, + co is 


2 


oT (53) 


Py, (2, )  2p,¢, aT os 
which for the metals considered is in the range of 10 to 
50 kilobars (see Table I). Certainly the equations of 
state derived above are not valid in the neighborhood of 

C (4.e., where Ch + oo) or in the region (see Fig. 2) to 
the left of C (1.e., where ¢50). Numerical computations 
for the 16 metals show that for pressures up to at least 

5 megabars (and in most cases much higher) the pressure 
increases monotonically with increasing z and decreasing V 
on the curve C. Furthermore for P<P,,(2, ) and P>P,,(z, ), Cc 
lies respectively to the right and to the left of the 
Hugoniot (as shown in Fig. 2). Along the Hugoniot, the 
specific heat CoH rapidly decreases from oo at Py(Z ) to 
a minimum (at which the pressure is at most 0.2 megabars 


as can be seen in Table II) and then slowly increases. 
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The anomalous behavior of Cy leads to the conclusion 
that it is impossible to postulate the simultan<-vous 
existance of (1) a Hugoniot satisfying Eq. (3) and (14) of 
isentropes satisfying the mirror-image assumption, both in 
the low pressure neighborhood of P.,V,,E,, and (441) of a 
constant thermal coefficient of volume expansion along 
P<O, Walsh, et al, using the Mie-Griineisen equation 
prescribed in addition to (1) and (1i1), the experimental 
values for cy along PeO. Since the Gruneisen ratio was 
taken as a function of V only it was easy to show that Cy 
was a function of S only. Therefore Cp was well-behaved 
in their calculations. The conclusion is that the mirror- 
image assumption in the neighborhood of the low pressure 
Hugoniot given by Eq. (3) leads to the anomalous behavior 
of Cp: However it is just in this region that it is not 
possible to determine experimentally a unique Hugoniot. 
For sufficiently low pressures the effects of material 
rigidity (which have been neglected in the above analysis) 
give rise to elastic-plastic wave structures. Here an 
arbitrary assumed shock wave will break up into two com- 
pression waves, an elastic wave which precedes a plastic 
wave. This "two-wave" structure indicates that at some point 
on the Hugoniot locus of P,V states, there is a violation of 
the Bethe-Weyl necessary condition (°P/av*),>0, for the 


existance of a stable shock in a single phase fluid. 


16 
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A. Special Results on the Hugoniot 


It is readily seen that on the Hugoniot curve, 


Pa = Poe z(z+1)/a, (54) 

Vy/Vy = 2 , z/{a(z+1)], (55) 

BE, = EB, + $ c6z*/a®, (56) 

T, = T, F(z)(1+2/(2+1)], (57) 

U = ¢,(z+1), (58) 
ac’z(z+1)7¢(z)/a 


“pH” z-eaaT (z+1)g(z)F(z)’ (59) 


c= e,(2z+1)(1-z/[a(zt1)]), (60) 


17 
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v. RESULTS FOR 16 METALS 


In Sections III and IV it was shown that under the 
mirror-image assumption a complete thermodynamic descrip- 
t4on of a metal depended only on the experimental values of 
Co and a, and the coefficient of volume expansion along the 


isobar, PaO. These data are listed in Table I. The U, Uy 


7 by apply- 


curves of the 16 metals of Table I were obtained 
ing the method described at the very beginning of Section I. 
Different known shock velocities in a brass standard, whose 
equation of state is of the Mie-Gruneisen form, give rise to 
different measured shock velocities in the specimen (any one 
of the 16 metals), and hence to a U,u, locus for the speci- 
men. Columns 7 and 8 of Table I compare the specific heat 

eC. = ac“/a at PosVy (see Eq. (47)) obtained from the mirror- 


p 


image assumption with the handbook values o (3) listed by 


Walsh, et al.? : 

A summary of calculations and results based on this 
assumption is given in Table II and these results may be 
compared with those obtained by McQueen and Marsh using the 
Mie-Griineisen equation of state (see Table IV of Ref. 4). 
It is found that the mirror-image assumption always leads 
to smaller values of the free-surface volume (1.e., at 


PxO) resulting from expansion from the shocked state. For 


example, for copper, corresponding to the entries of column 


18 
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10 of Table II, the Mie-Griineisen equation gives for the 
free-surface volumes, WV, = 1.000, 1,.000,. 1.001, 1,005, 
1,006, -14010,. 1<015,. 1,019, 1,025; 2.051, 1.057; 1.044, 
1.051, 1.058, 1.066, and 1.079. Furthermore, along the 
isentropes tabulated, WV, is always smaller for identical 
pressures except in the case of thallium (see columns 14 
and 15 of Table II); and for identical P the differences 
between the compressions, (V(P,z)-V,(2)]/Vy(z), calculated 
with each equation of state (starting from the same Pu(z)s 
Viy(Z))s are smallest at high pressures (particularly for 
Mo, Ni, Ti, V, and W) and greatest at low pressures. As 

a consequence the temperatures, both along the Hugoniot and 
along the isentropes, are also lower (and for some metals 
very much lower) with the mirror-image assumption, For 
example, for copper, corresponding to the entries of 
colum 4 of Table II, the Mie-Gruneisen equation gives for 
the Hugoniot temperatures (in °K), Ty + 293, 334, 391, 472, 
582, 719, 881, 1068, 1277, 1506, 1755, 2020, 2301, 2596, 
2902, and 3042. 

The results of these computations demonstrate that 
the thermodynamic description obtained by the mirror-image 
assumption causes the metal to appear less compressible 
than when the thermodynamic description is obtained from 


the Mie-Griineisen equation of state, 
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Table I. Input data and comparison with zero-pressure data. The values of 
Por Sor a and of) are taken from Reference 4 (see their Table III); and the thermal 
coefficient of volume expansion a@ and the specific heat (3) are handbook values from 
Reference 3 (see their Table I), except as indicated. The quantity e,(0,0) = oO ce/a 
appearing in  colum & is the theoretical specific heat obtained by letting 
z~0 in Eq. (47). The quantities o(B), of), and op (0,0) are values at the initial 
state P.w0, V). Pilz, )s the calculated shock pressure at which c.~ + 00, 18 conlputed 


from Eqs. (51) and (54), 


P 


bods, MEA he ee o(®) (3) 0,(0,0)  Pulzy) 

Material g/cm? /°K cm/sec om/usea (cal/g)/K (ca1/e)°x Kiloba: 
Ag 10.49 56.7 1.586 0.3243 0.319 0,056 0.090 40.9 
Au 19.24 42.6 1.560 0.3075 0.305 0.031 0.043 49.2 
ca 8.64 89.4 1.671 0.2443 0.241 0.055 0.076 32.6 
Co 8.82 36.9 1.330 0.4748 0.463 0.099 0.149 45.6 
Cr 7.10 18.6 1.465 0.5217 0.515 0.111% 0.083 21.7 
Cu 8.90 49.5 1.497 0.3958 0.398 0,092 0.124 4h 2 
Mo 10.20 15 1.238 0.5157 0.519 0.061 0.077 2k 4 
Ni 8.86 39 1.445 O.4646 0.463 0.105 0.139 46.8 
Pb 11.34 85.1 1.517 0.2028 0,202 0.030, 0.055 27.3 
Sn 7.28 60 1.476 0.2640 0.276 0.054 0.068 19.9 
Th 11.68 33.3 1.278 0.2132 0.205 0.030 0,028 10.9 
T4 4.51 25.5 1.089 0.4779 0.484 0.126 0.128 15.9 
Tl 11.84 124 1.515 0.1859 0.183 0.031 0.062 34.2 
6.1 23.45 1,210 0.5108 0.518 0.115% 0.121 22.6 

19.17 17.45 1.268 0.4005 0.405 0.032 0.053 32.2 

Zn 7.14 100 1.559 0.3050 0,303 0.092 0.143 47.5 


S Smithsonian Physical Tables, prepared by W. E. Forsythe (Smithsonian Institution, 
“Washington, D.C., Ninth Revised Edition, 1954), pp. 155, 157 
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Fig. 1. Pressure versus particle velocity curves 
for an initial shock in a metal standard followed by a 
reflected shock or rarefaction wave, and a typical graphi- 


cal solution to determine PysUy for a test specimen. 
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Fig. 2. Pressure versus volume curves for a metal 
that satisfies the mirror-image assumption. The dash«dot 
curve, C, is defined by Eq. (50) and is the locus of 
points on which Cp" + © depending upon the direction of 
approach. The intersection of C with Hugoniot occurs at 
Pi(Z )s Viz (zy ) where Py(z, ) lies between 10 to 50 kilobars 
for the 16 metals considered. The figure is split into a 


megabar region and a "less than 60 kilobar region", 
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